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Abstract. We consider the propagation of polarized photons in optical fibers under the action 
of randomly generated noise. In such situation, the change in time of the photon polarization can 
be described by a quantum dynamical semigroup. We show that the hierarchy among the decay 
constants of the polarization density matrix elements as prescribed by complete positivity can be 
experimentally probed using standard laboratory set-ups. 



1. Introduction 

A quantum system immersed in a large environment represents the paradigm 
of a physical situation that is often encountered in quantum optics, quantum 
chemistry and atomic physics. Since it is not isolated, such system is called open: 
it exchanges energy and entropy with the external environment, whose dynamics is 
typically unaffected by the presence of the small subsystem. The evolution of such 
open quantum systems, obtained by tracing over the (infinite) degrees of freedom 
of the environment, is in general very complicated, exhibiting non linearities and 
memory effects. 

However, when the coupling between the two systems is weak, a very com- 
mon situation in many practical applications, the reduced time evolution can be 
described by a one-parameter (=time) family of linear maps {7*}, a so-called quan- 
tum dynamical semigroup PP-[H]- These linear transformations map states of the 
subsystem, described by density matrices, into states, while possessing many desir- 
able physical properties. In particular, these maps exhibit irreversibility, encoded 
in the semigroup property (i.e. in the forward in time composition law), decoher- 
ence and dissipation, the typical noisy effects induced by the environment. Further, 
these generalized time evolutions turn out to be completely positive, a property 
that guarantees the physical consistency of the dynamics in all physical situations. 

Nevertheless, in many phenomenological treatments of the dynamics of open 
quantum systems, this last property is often dismissed as irrelevant, nothing more 
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than a mathematical artifact [H]-jlL)j. The argument supporting this point of view 
is based on the following considerations. In order to represent a physical state 
of the subsystem S, a density matrix p must be a positive operator, since its 
eigenvalues have the meaning of probabilities; this is at the root of the statistical 
interpretation of quantum mechanics. The time evolution ^(O) i— ► pit) = 7t[p(0)] 
must then preserve this fundamental property, and therefore map a positive initial 
p(0) into a positive final pit). Such a property of the linear transformation jt is 
called positivity; it is apparently sufficient to assure the physical consistency of the 
dynamics. 

On the other hand, complete positivity is a more restrictive requirement than 
simple positivity: it guarantees the positivity not only of jt but also of the dy- 
namics of a larger system built with the system S statistically coupled to a second 
ancillary system S', which however remains inert under the time evolution. The 
dynamics of this enlarged system is then described by 7j = 7$ <8> id, where id rep- 
resents the identity operation. Positivity of Ij for any ancillary system S' means 
complete positivity of the map jf It is worth noting that this property is inti- 
mately related to entanglement, i.e. to the possibility that S may have interacted 
in the past with S' and become quantum correlated with it [H lllj. 

Admittedly, although logically stringent, the above motivation of why an open 
system dynamics should be completely positive is not very appealing from the 
physical point of view; indeed, the ancilla system S' is remote from S and fixed. 
Nevertheless, a more concrete scenario can be offered, that of two equal, mutually 
non interacting systems S, immersed in the same external bath; such a situation 
is not uncommon and it is for instance encountered in high energy particle physics 
in the study of correlated neutral mesons |12| 1131 114j . In this case the two systems 
evolve with the product map A = 7t <8> 7*j and one can show that the positivity of 
A is equivalent to the complete positivity of jt 

Another important advantage of having a completely positive open system 
dynamics is that this property fully characterizes the form of the map |16M17| . In 
particular, when S is a two-level systems, as discussed below, the decay times of the 
diagonal (A) and off-diagonal (T2) elements of the corresponding density matrix 
are seen to satisfy a characteristic order relation: 2 A > A Pd]. It is precisely 
the presence of such hierarchy that is often questioned in the phenomenological 
literature on open quantum systems. 

The aim of the present investigation is to analyze in detail the possibility of 
a direct test of the above inequality in the laboratory, using polarized photons. 
The basic idea is to study the change in their polarization while the photons travel 
across a weakly coupled, noisy environment. A possible practical realization of this 
scenario involve photons with a given polarization injected inside a high quality, 
polarization preserving optical fiber, which is subjected to random noise generated 
by external high frequency sound waves; as we shall see, the random noise will play 
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a role similar of that of an external environment. By making the photons traverse 
the fiber twice with the insertion of a Faraday mirror at one end, it is possible 
to measure the ratio between the diagonal and off-diagonal elements of the 2x2 
density matrix representing the final photon polarization state, and thus determine 
the combination 2T\ — T<i- 

All the necessary techniques needed for the realization of the apparatus and the 
realization of the measure are standard and available in modern quantum optical 
laboratories. We hope that our results will stimulate the actual realization of the 
experiment, thus clarifying in a controlled setting the role of complete positivity 
in open quantum dynamics. 



2. Master equation 

In describing the evolution of polarized photons in a optical fiber we shall 
adopt the standard effective description in terms of a two-dimensional Hilbert 
space, the space of helicity states Any vector in this space represents a 

given polarization and can be identified by two angles 9 and (p: 

\6,<p) = cos 6\+) + e ilp sm9 |-) , (1) 

where |+) and |— ) are two orthonormal basis vectors, representing linearly po- 
larized states. Another convenient basis in this space is given by the circularly 
polarized states: 

\R) = -A=(|+) + , \L) = -4(1+) - • (2) 



With respect to this basis, any (partially) polarized photon state can be rep- 
resented by a 2 x 2 density matrix S; as already pointed out before, this is a 
hermitian, positive operator, i.e. with positive eigenvalues, normalized to have 
unit trace. 

The time evolution of the photons inside the fiber, while subjected to random 
noise, can be cast in a standard Liouville - von Neumann form: 

^jp- = -i [Ho ,m} + L t m)} ■ (3) 

The first piece on the r.h.s. describes the propagation of the photons in absence of 
noise; in the chosen basis, the effective hamiltonian Hq can be cast in the general 
form: 

H = fn-a, (4) 

where a = (a%, 02, 03) is the vector of Pauli matrices. For sake of generality, we 
have kept ujq nonvanishing, in order to take into account possible birefringence 
effects due to the fiber mechanical bending. 
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The fiber is further subjected to random noise, generated by sound waves of 
frequency much higher than the inverse flight time of the photons in the fiber. As 
a consequence, the travelling photon will propagate in wave guide that is randomly 
changing; it is as if it was effectively moving into a rapidly fluctuating medium, 
which can be described by classical stochastic fields. The action of this fields on 
the travelling photons can then be expressed via the commutator with a time- 
dependent hermitian matrix F(t), 

Lt\E(t)] = -i[F(t),V(t)], F(t) = F(t) ■ a , (5) 

whose components Fi(t), F2(t), Fs(t) form a real, stationary Gaussian stochas- 
tic field F(t); in general, they have nonzero constant means and translationally 
invariant correlations: 

Gij(t -s) = (F^Fjis)) - (Fi(t)) (Fj(8)) , i,j = 1,2,3 . (6) 

Since the generalized hamiltonian F(t) in © involves stochastic variables, the 
density matrix solution of the equation of motion (j2J), is also stochastic. 

Instead, we are interested in the behaviour of the reduced density matrix p(t) = 
which is obtained by averaging over the noise; it is p(t) that describes the 
effective evolution of the photons in the randomly behaving fiber. By making the 
additional assumption that photons and noise be decoupled at t = 0, so that the 
initial state is p(0) = (S(0)) = S(0), a condition very well satisfied in practice, 
an effective master equation for p(t) can be derived by going to the interaction 
representation, where we set: 

= e ltH ° S(t) e- itHo , a(t) = e itH ° a e~ itH ° , L t [ ] = -i[P{t) ■ a{t), ' . 

(7) 

The time evolution of the reduced density matrix in this representation, p(t) = 
can then be expressed as a series expansion involving multiple correlations 
of the operator Lt. 

p(t)=M t [p(0)] = jr td Sl fW" r^dsk (L Sl L S2 ---L Sk )[p(0) \ ■ (8) 
jZr Q J0 Jo Jo 

The sum Mt[p] = X)a^=o ^k\P] i n © with Nq[p] = p can be formally inverted, and 
a suitable resummation gives (a dot represents time derivative) |24| : 



dp(t) 
dt 



M t Mf l [p(t)] = {iVi + (N 2 -N 1 N 1 ) + .. .}{p(t)] . (9) 



Since the interaction of the travelling photons with the stochastic field is weak, 
one can focus the attention on the dominant terms of the previous expansion, 
neglecting all contributions higher than the second-order ones. Further, since the 
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where 
and 



characteristic fluctuation time of the random noise is by assumption much smaller 
than the typical photon time of flight across the fiber, the memory effects implicit in 
© should not be physically relevant and the use of the Markovian approximation 
justified. This is implemented in practice by extending to infinity the upper limit 
of the integrals appearing in and iVW H|-|3|. 

By returning to the Schrodinger representation, one finally obtains 

d -P^ = -i[H,p{t)] + L[p(t)] , (10) 
H = H + H x + H 2 = Co ■ a , (11) 

1 3 

L \p\ = 2 J2 Ci i [ 2a iP a i ~ {°i a 3 > P}\ ■ ( 12 ) 
*,i=i 

The effective hamiltonian H differs from the one in absence of noise Hq by 
first order terms (coming from the piece in Q) depending on the noise mean 
values: 

H x = (F(t)) ■ a , (13) 
and by second-order contributions (coming from the second-order terms in Q), 

3 

H2 = tijk Cijak , (14) 

i,j,k=l 

involving the noise correlations (jHJ) through the time-independent combinations: 

3 poo 

dj = Y, d t Gik{t) U kj (-t) , (15) 
k=i Jo 

where the 3x3 orthogonal matrix U(t) is defined by the following transformation 
rule: e ltH ° Oi e~ ltH ° = J2j=i U%j(t) aj. On the other hand, the contribution L[p] 
in (|12j) is a time- independent, trace-preserving linear map involving the symmetric 
coefficient matrix Cjj = Cjj + Cji, which needs to be positive, Cjj > 0, in order 
to generate a completely positive dynamics Q]-|l]- It introduces irreversibility, 
inducing in general dissipation and loss of quantum coherence. 

Altogether, equation (|1()| ) -(|12| ) generates a semigroup of linear maps, r t : 
p(0) 1 — ^ p(t) = lf[p(0)], for which composition is defined only forward in time: 
i~i o r s = with t, s > 0; it is usually referred to as a quantum dynamical 

semigroup 
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It is convenient to expand the photon density matrix in terms of the Pauli 
matrices and the identity ctq, 



1 

P= 2 



a + p-a 



(16) 



Then, the linear equation Q1U|) reduces to a diffusion equation for the components 
Pi > Pii P3 of the (Bloch) vector p: 



dm 

dt 



-2Hp(t) ; 



(17) 



the entries of the 3x3 matrix 7i can be expressed in terms of the coefficients uJi 
and Cij appearing in the hamiltonian and noise contribution in Ullj) , ()12j) ^ |1] : 



U 



a 

b - w 3 

C + Ll>2 



b + iv 3 
a 



C-0l> 2 
P + Ul 
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(18) 



with a = C 22 + C 33 , a = C n + C 33 , 7 = Cn + C 22 , 6 = -Ci 2 , c = -C i3 , /3 = -C 23 . 
The condition of complete positivity, Cij > 0, can then be expressed more explicitly 
through the following inequalities: 



2R = a + 7 
25 = a + 7- 
2T = a + a - 



-a > 
a > 
•7 > 



RS>b 2 , 
RT>c 2 , 
ST>(3 2 , 



RST >2bc(i + R(3 2 + Sc 2 + Tb 2 . 
The solution of (|17|) involves the exponentiation of the matrix 

-2Ht 



p(t) = M(t) p(0) 



M(t) 



(19) 



(20) 



From this relation, one immediately sees that the 3x3 matrix M(t) represents 
the (reduced) Mueller matrix connecting the initial Stokes vector p(0) with the 
evolved one p(t) at time t [TH l22| l23 | l^fi] . 



3. Complete positivity 

As shown in the previous section, the entries of the matrix 7i in ()18|) are directly 
related to the stochastic field correlations in Q. For a generic noise, one expects 
all parameters in H to be nonvanishing; only physical considerations may allow a 
simplification and therefore a manageable, explicit expression for p(t). 

For a typical random noise, one can generically assume the correlation functions 
to have an exponentially damped form. Further, the off-diagonal correlations are 
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usually much suppressed with respect to the diagonal ones, so that, without much 
loss of generality, one can write: 

G ij (t) = G i e- Xi 1*1 Sij , G u Xi>0. (21) 

Using the expression (|3J) for the hamiltonian Ho, one easily finds the form of the 
matrix Uij(t) giving the time evolution of the Pauli matrices, 

3 

Uij(t) = TiiUj + (5ij — riirij) cos uiot — ^ £ijk n k sincjo^ (22) 

k=i 

and then through the definition Q15|) . one finally obtains: 



Cij — AjAj 



2 , , 3 

(23) 



r ^0 , ^0 \ - 



fe=l 



where 

A, = ^f-J, • (24) 
A t + W 

The symmetric and antisymmetric pieces of Cij give the noise contributions to the 
dissipative ()12j) and hamiltonian (|14|) parts of the evolution equation (|10jl . 

In order to further simplify the treatment, we shall assume the unit vector n 
defining the starting hamiltonian Hq to be directed along the third axis and the 
stochastic field F(t) to have zero mean. In this case, the effective hamiltonian in 
(fllft remains proportional to 03 , H = ua3 , with a new frequency u> = 1V3 explicitly 
containing zero and second order pieces: uj = u>o/2 + uq (Ai + A2). Similarly, also 
the dissipative contributions simplifies, since one finds c = (3 = 0, while 

a = 2A 2 A 2 + 2^, 6 = wo(A 2 -Ai), 
A3 

a = 2AiA! + 2^, 7 = 2A1A1 + 2A 2 A 2 . (25) 
A3 



Notice that the conditions assuring the complete positivity of the time evolution 
given by the inequalities in (|19[) are not all automatically satisfied by the as- 
signments in (|25jl. The following relation needs to be imposed, Wq(^-2 — Ai) 2 < 
4A1A2A1A2: it can always be fulfilled by a suitable choice of the noise parameters 
in ((UJ. 

This situation is not exceptional: the derivation of a physically consistent, 
Markovian reduced dynamics starting from the exact Liouville - von Neumann 
equation in (j3J) is in general very involved and can be treated with the necessary 
mathematical rigor only in special cases ^-(Sl- In general, as briefly discussed 
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in Section 2, one instead resorts to various approximations, justified by physical 
considerations. However, this naive procedure does not always lead to physically 
consistent reduced dynamics [3j: further conditions, like that of complete posi- 
tivity, need to be imposed at the end. This explains why this property is often 
dismissed as irrelevant in phenomenological applications, and therefore why it is 
so important to verify its fulfillment in an experimentally controlled setting. 

With the above assignments, the exponentiation of 7i can be easily evaluated 
in closed form, and for the Mueller matrix one explicitly finds: 



M(t) 



e -(a+a)t 
e -(a+a)t B _fy 



e -(a+a)t B ^ 
e -(a+a)t a_ 








-27* 



(26) 



where 



ct — (X 

A±(t) = cos(2flt) ± — — sin(2Qi) 



B±(t) 



b±uj 

n 



20 
sin(2fii) , 



(27) 



with O = [uj 2 - b 2 - (a — a) 2 /4] 1//2 . The evolution matrix in 1)26(1 contains both 
oscillating and damping terms, and can be considered as a generalization of a 
rotator |21|-|23|. Notice, however, that the oscillator behaviour depends on the 
magnitude of effective frequency ui with respect to the dissipative parameters a, b 
and q; when u < [b 2 + (a — a) 2 /4] 1//2 , the frequency Q becomes purely imaginary 
and M{t) contains only exponential terms. 

In any case, for large times, i.e. for infinitely long fibers, the dumping terms 
always dominate |27j and the Mueller matrix becomes that of a total depolarizer; 
in other words, as result of the action of the stochastic noise, the photon den- 
sity matrix p(t) will become asymptotically proportional to the unit matrix (i.e. 
p(t) ~ 0), independently from the initial state p(0). 

Through the definition (|16|) . the evolution matrix (|26|) gives the time be- 
haviour of the entries of the photon density matrix p(t), which in turn can be 
experimentally determined using suitable tomographic procedures. Therefore, at 
least in principle, thanks to the different time-dependence of the entries in (|26() . 
one can measure the magnitude of all the dissipative parameters a, b, a and 7 
and thus check the conditions of complete positivity, that in particular requires: 
a + q — 7 > 0. 

Unfortunately, this can not be done by a single passage of the photons in 
the noisy fiber, since in this case the time t appearing in (|26|) is fixed, being the 
time of flight along the fiber. In order to isolate the exponential damping terms 
from the oscillating ones, t needs to be varied. This can be achieved by placing 
at the end of the fiber a Faraday mirror [2^, which inverts the polarization of 
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the photons, while reflecting them back into the fiber. In this way, travelling 
backward, the rotation induced by the standard hamiltonian Hq on the photon 
polarization is "undone", while the dissipative effects due to the stochastic noise 
further accumulates. Indeed, the backward evolution is given by a Mueller matrix 
«M(i) still of the form (|26|). but with the frequency ujq replaced by — wo, which in 
turn gives to — > —uj and b — > —b, while leaving the remaining dissipative parameters 
unchanged. 

Thanks to the semigroup property of the dynamics, the complete evolution 
for the double passage of the photons in the fiber is obtained by composing the 
two Mueller matrices, and therefore, after having travelled for a time It inside 
the noisy fiber, the Stokes parameters representing the photon polarization state 
become at the end: 

p(2t) = M{t) -M{t)p(0) . (28) 

Let us now consider an initially linearly polarized photon, so that p(0) = 
p( + >(0) = (1,0,0). After having passed once trough the noisy fiber, this Stokes 
vector becomes p^ + \t), as in (|2"U)) . while further evolves to p( + \2t), as in (|2"5|) . 
after having travelled backwards to the beginning of the fiber. Similar results 
hold for a circularly polarized initial photon, for which: p(0) = p( R '(Q) = (0,0, 1). 
One can then form the following combination of components of the above Stokes 
vectors: 

I A Kt). 



n(t) 



A%) 



(29) 



Using the explicit expressions for these Stokes components as obtained from (|26|) . 
one easily finds: TZ(t) = exp[— 2(o + a — 7) t]. Recall that complete positivity 
requires the combination a + a — 7 to be nonnegative. Therefore, by measuring 
the components of the Stokes vectors appearing in (|29fl . one can determine the 
quantity lZ{t) and thus check whether or not 1 

n(t) < 1 . (30) 

Although this condition is in general only necessary for complete positivity, it 
becomes also sufficient for the phenomenologically relevant situation for which 
Ai = A2. Recalling (|25|) . in this case one has 6=0 and further a = a, so that 
2a > 7 is the only surviving inequality of those listed in ()19|) . 

The two nonvanishing dissipative parameters a and 7 have now the meaning 
of inverse relaxation times for the off-diagonal and diagonal entries of the photon 
density matrix p{t) pQ. They are usually called I/T2 and 1/Ti, respectively. The 
condition that assures the complete positivity of the open dynamics for photons 

1 Note that this inequality determines the sign of the combination a + a — 7 for any t, i.e. 
irrespective of the length of the fiber, provided it behaves in a stochastic way. This fact might 
help reduce the systematic uncertainties in the actual experimental test of 1301 . 
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travelling along the noisy fiber is therefore precisely 2T\ > T2, as was mentioned 
in the introductory remarks. 

In summary, we have shown that the condition assuring the complete positivity 
of the dynamics of polarized photons in a noisy fiber can be cast in the form (|3U|) , 
By measuring the combination 1Z in (|29|) . it can be tested using set-ups and tech- 
niques that are routinely used in optical laboratories. We find this possibility very 
intriguing and hope will trigger the interest of the vast community of specialists 
working in quantum optics experiments. 
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